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Abstract We briefly review recent developments in 
black hole accretion disk theory, emphasizing the vital 
role played by magnetohydrodynamic (MHD) stresses 
in transporting angular momentum. The apparent uni- 
versality of accretion-related outflow phenomena is a 
strong indicator that large-scale MHD torques facilitate 
vertical transport of angular momentum. This leads to 
an enhanced overall rate of angular momentum trans- 
port and allows accretion of matter to proceed at an 
interesting rate. Furthermore, we argue that when ver- 
tical transport is important, the radial structure of the 
accretion disk is modified at small radii and this af- 
fects the disk emission spectrum. We present a sim- 
ple model demonstrating how energetic, magnetically- 
driven outflows modify the emergent disk emission spec- 
trum with respect to that predicted by standard accre- 
tion disk theory. A comparison of the predicted spec- 
tra against observations of quasar spectral energy dis- 
tributions suggests that mass accretion rates inferred 
using the standard disk model may be severely under- 
estimated. 

Keywords accretion, accretion disks — black hole 
physics — (magnetohydrodynamics:) MHD — radia- 
tion mechanisms: general — turbulence waves — galax- 
ies: active — X-rays: binaries 



1 Introduction 

The standard theory of astrophys ical disk accretion was 
formulated over thirty years ago ( Pringle fc Reed 1972 : 



Zdenka Kuncic 

School of Physics, University of Sydney, Sydney, NSW, Australia 
2006 

Geoffrey V. Bicknell 

Research School of Astronomy &; Astrophysics, Australian Na- 
tional University, Cotter Rd., Weston ACT, Australia 2611 



Novikov fc Thornd Il973t IShakura fc SunvaevI llQTSi ). 

Since then, arguably the most important theoretical 
milestone has been the demonstration via numerical 
simulations that the removal of angular momentum re- 
quired for accretion to proceed is mediated by magneto- 
hydrodynamic (MHD) turbulence driven by th e weak- 
field magnetorotational instability (MRI - see Balbuj 
20031 for a review). Notwithstanding these ground- 



breaking developments, it still remains unclear precisely 
how energy can be channelled from the bulk accretion 
flow to a diffuse region outside the disk where high- 
energy processes can operate. This outstanding issue 
is inextricably linked with the formation and ubiquity 
of outflow phenomena associated with accretion disks. 
Indeed, the formation and nature of relativistic jets re- 
mains one of the most formidable theoretical problems 
in this fleld. 

In this paper, we highlight some recent theoretical 
progress made towards understanding vertical t rans- 
port in accretion disks (jKuncic fc Bicknell 120041) . A 
macroscopic, mean field approach is adopted for the 
MHD turbulence. This allows us to identify the dom- 
inant mechanisms responsible for angular momentum 
transport as well as the main contributions to the 
global energy budget of the system. We also demon- 
strate quantitatively how outflows, both Poynting-flux- 
dominated and mass-flux-dominated, can modify the 
spectrum of disk emission substantially from that pre- 
dicted by standard accretion disk theory. 



2 MHD disk accretion 

Here, we present statistically-averaged equations in a 
cylindrical (r, 0, z) coordinate system for a fluid which 
is time-independent and axisymmetric in the mean 
{{d/dt) ~ (d/dcj)) ~ 0, where the angle brackets de- 
note ensemble-averaged quantities). We use Newto- 
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nian physics throughout, with the gravitational poten- 
tial (/)g = ~GM{r^ + z2)-i/2, xhc fluid description 
is valid down to the innermost marginally stable orbit 



Gr^, where 



GM/c is the gravitational radius 



of a black hole of mass M. The mean-field conserva- 
tion equations are integrated vertically over an arbi- 
trary disk scaleheight, h = h{r). Quantities calculated 
at the disk surface (z = ±h) are denoted by a ± super- 
script, , and we assume refiection symmetry about 
the disk midplane, so that l^"*"! = \X~\. Midplane 
values of variables are denoted by Xq. We assume h 
is much less than the radius (the "thin disk" approxi- 
mation), so that quantities of order h/r and dh/Ar are 
neglected. 

We adopt a mass-weighted statistical averaging ap- 
proach in which all variables are decomposed into mean 
and fluctuating parts, with intensive variables such as 
the velocity mass averaged according to 



Vi = U, + Vi 



{pv[) = 



(1) 



while extensive variables, such as density, pressure and 
magnetic fleld, averaged the following way: 



P = P + P 
p — p + p' 
Bj = Bj + B'i 



ip') - 

ip') = 

(Bl) = 



(2) 
(3) 
(4) 



Note that intensive averages are denoted by a tilde, 
extensive averages are denoted by a bar and fluctu- 
ating components are denoted by a prime. The fluc- 
tuating velocity components are restricted to subsonic 
speeds because the MRI is a weak-field instability which 
drives subsonic turbulence with (pv'^) < (pc^) ^ pv^ ~ 
pr^l^K^, where Cg = {kT/ pmpY^^ is the local sound 
speed and f^K = {GM/r^y^^ is the keplerian angular 
velocity. The only restriction we place on the mean fluid 
velocity components is that they satisfy "D^ 3> Vr,Vz 
and that Vr and are independent of z, simplifying 
the vertical integration. The combined magnetic and 
Reynolds stresses are defined by 



(5) 



In what follows, we summarize the salient equa- 
tions for radial and vertical transport of mass, an- 
gular momentum and energy. A full derivation of 
these equations can be found in our earlier paper 
( Kuncic fc Bicknelllliool . For clarity, we have simpli- 
fied the presentation of the relevant equations by omit- 
ting negligible correlation terms in the energy equatior0 



^In particular, triple correlation terms of the form {Ujv'j) are 
negligible compared to analogous correlations with the mean fluid 
velocity {tij)vj. 



and by removing the notation for averaged extensive 
and intensive quantities. Thus, averaged quantities are 
implicitly assumed. 

2.1 Mass transfer 

Vertical integration of the mean-field continuity equa- 
tion gives 



— / 2iTrpVr dz -\- Anrp'^v^ = 
d?- J-h 



(6) 



We now introduce the usual definitions for the surface 
mass density, 



r+h 

I](r) = / pdz 
J-h 



(7) 



and mass accretion rate, 

Afa(r) = 27rrS(-v,) . (8) 
We also introduce an analogous mass outflow rate. 



(9) 



associated with a mean vertical velocity at the disk 
surface, i.e. at the base of a disk wind. In terms of the 
above definitions, the vertically integrated continuity 
equation, ([5]), can be written as 



dr dr 
implying that 



(10) 



Afa(r) + M^{r) = IvU{r{) + Afw(n) = constant = M 

(11) 

where Afa(n) + M^{r\) is the total mass fiux at the 
innermost stable orbit, n. Equation (fTO|) impHes that 
under steady-state conditions, the radial mass inflow 
decreases towards small r at the same rate as the ver- 
tical mass outflow increases in order to maintain a con- 
stant nett mass flux, M , which is the nett accretion 
rate at r = oo, i.e. M = M^{oo). 

2.2 Angular momentum 

The azimuthal component of the momentum equation 



1 ^ ^ 2 \ / 
^— (r pvrv^) -t- — {pv^v, 

1 d 



(12) 



p2 Qj. 



d{t4,z) 
dz 
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Integrating this equation over z and using the mass 
continuity equation (|10p gives 



d 
dr 



where 



Tr4, — I {tr<ji) 



(13) 



(14) 



is the vertically integrated r0 stress. The terms on the 
the left hand side of eqn. describe radial trans- 
port of angular momentum associated with radial in- 
flow (accretion) and MHD stresses acting over the disk 
height. The terms on the right hand side describe verti- 
cal transport of angular momentum resulting from mass 
loss in a wind and MHD stresses on the disk surface. 
The magnetic part of the MHD stresses are given by 
{tr4,) - {BrB^)/47r and - {B^B,)/4:Tt (cf eqn.p 

Radially integrating the angular momentum equa- 
tion (fO)) gives 

Ma,v^r - Mi,{ri)v^{ri)ri = + 2TTrfTr4,{n) 



dr 



47rr 



dr . (15) 



This is a generalized conservation equation for angular 
momentum in accretion disks. It can be equivalently 
expressed in terms of the flux of angular momentum, 
J, as follows: 



L{r) - ja(n) jr{r) - Jrin) + (r) 



(16) 



where Ja = M^v^r is the angular momentum flux of 
the accreting matter, = — 27rr^Tr0 is the radial flux 
of angular momentum carried by the rcj) stresses and 
Jz — J^. ■yi^r dMa/dr — 47rr^(i02)+ dr is the vertical 
flux of angular momentum carried by outflowing matter 
and by the (f>z stresses. Thus, both radial and vertical 
transport contribute to the nett rate at which matter 
loses angular momentum and moves radially inwards. 
In other words, the mass accretion rate. Ma, is deter- 
mined by the nett rate at which angular momentum is 
transported radially outwards by MHD stresses acting 
over the disk height and vertically outwards by a mass 
outflow as well as by MHD stresses acting over the disk 
surface. 

The most luminous accretion-powered astrophysi- 
cal sources are inferred to be accreting matter at 
very high rates. Quasars, for instance, appear to be 
accreting close to or even exceeding the Eddington 



rate, Msdd = AirGMm^/ {-qa^c) « O.3 77o;i[Af7M0 yr'^, 
where rj — O.lryo.i is the efficiency with which gravi- 
tational binding energy is converted to radiation and 
M = \{f MjA4q is the mass of the central black hole. 
Identifying the dominant mechanisms responsible for 
such high mass accretion rates is a major challenge con- 
fronting accretion disk theory. 

It is interesting to compare the rate of angular mo- 
mentum transport indicated by ()15|) with that pre- 

diete d by standard accretion disk theory ( Shakura fc Sunvaev^ 
19731 ). In a standard disk, angular momentum is trans- 
ported radially outwards only and the stresses respon- 
sible for this are parameterized by a dimensionless pa- 
rameter a such that 



Trrb — acMYiT 



dr 



(17) 



In addition, the stresses are assumed to vanish at ri and 
Afa is necessarily constant with r (since there is no mass 
outflow). Thus, the angular momentum conservation 
equation, p^ . reduces to 



dr 



(18) 



The parameter a lies in the range < a < 1 for stresses 
resulting from subsonic turbulence. Values of a ap- 
proaching unity are required to account for the high 
accretion rates inferred in the most powerful sources, 
such as quasars. However, because the standard a— disk 
model remains a phenomenological prescription, it is 
unclear whether values a 1 correspond to a physically 
plausible realization of the actual turbulent stresses in 
real accretion disks. 

Numerical simulations may be uniquely capable of 
addressing this issue. Simulations of accretion disks 
have demonstrated conclusively that turbulent MHD 
stresses can indeed remove angular momentum from 
matter, thus facilitating the accretion process (see 
BalbusI I2OO3I for a review). However, 3D simula- 
tions of MRI-generated MHD turbulence in accretion 
disks have so fa r been unable to produce high mass 
accretion rate s ( Hawlev 20001 Stone fc Pringlel 2001; 
Hawlev ev all boOlt iHawlev fc BalbusI I2OO2I ). Typi- 
cal values of a are a few x 10~^. On the other 
hand, high accretion rates are recovered from 3D MHD 
simulations when a large -scale, open mean magnetic 



field is explicitly iriclude d (jSteinacker fc Henningll2001 
Kigure fc Shibatal 120051) . These simulations and oth 



ers fe.g. ISalmeron et al.ll2007l) . as well as semi-analytic 
models (e.g. Campbell 2003f) . show that at small disk 
radii, vertical transport of angular momentum by a 
large-scale magnetic torque is more efficient than ra- 
dial transport by MHD turbulence. 
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The numerical results indicate that MHD turbu- 
lence, whilst important for the microphysics of accre- 
tion disks, cannot be solely responsible for the removal 
of excess angular momentum from accreting matter. As 
eqn. (|15p shows, angular momentum can also be re- 
moved by vertical components in the mean fluid and 
magnetic fields (i.e. nonzero Vz an d (B^B^,) compo- 
nents at the disk surface - see also iKonigl fc Pudritz 
Large-scale MHD effects in the form of mass out- 
flows and magnetic torques must therefore be primar- 
ily responsible for enhanced angular momentum trans- 
port in accretion disks. The fundamental relationship 
between angular momentum and energy (dE = QdJ) 
then also implies that large-scale MHD effects must also 
be responsible for the high-energy phenomena associ- 
ated with accreting sources which must necessary orig- 
inate outside the disk in a relativistic jet or magnetized 
corona. 

In the following section, we examine the global en- 
ergy budget of accretion disks in which angular momen- 
tum transport is prescribed by the generalized angular 
momentum conservation equation (jlSp . 

2.3 Radiative disk flux 

We now consider energy conservation in MHD disk ac- 
cretion. Again, we refer t he reader to our earlier paper 
( Kuncic fc Bickneli 2004 ) for details of the derivations. 

Accretion power is the rate at which gravitational 
binding energy is extracted from the accreting matter. 
This energy can be converted into mechanical (e.g. ki- 
netic, Poynting flux) and non-mechanical (e.g. radia- 
tive) forms. The rate at which this occurs is determined 
by keplerian shear in the bulk flow, Sr4, = ^rd^l/dr, 
with dn/dr = -|f^/r. The rate per unit disk surface 
area at which energy is emitted in the form of electro- 
magnetic radiation is determined by the internal energy 
equation. If there are negligible changes in the internal 
energy and enthalpy of the gas, then the disk radiative 
flux, , is approximately equal to the rate of stochas- 
tic viscous dissipation of the turbulent energy, which 
occurs on the smallest scales at the end of a turbulent 
cascade. If there is negligible transport of turbulent 
energy from the source region, then turbulent energy is 
locally dissipated at a rate equivale nt to its production 
rate, {tij)sij ~ {tr4,)sr4, (see § 2.5 in lKuncic fc Bickneli 
2OO4I for a discussion on the relative importance of pro- 
duction, transport and dissipation of turbulent energy) . 
The internal energy equation then implies 



the disk surface: 



m 

dz 



{tr 



Vertically integrating over the disk height yields the fol- 
lowing expression for the radiative flux emerging from 



dr 



(20) 



The level of the turbulent MHD stresses T^^ avail- 
able for internal dissipation depends on how efficient 
other processes are at converting the extracted accre- 
tion energy into other (mechanical and non-mechanical) 
forms. In other words, T^^ is specified by the angular 
momentum conservation relation ()15|) . which gives 



27rr2 
2 



1 - 



(7) Tr,{n) 



2717''- 



dr 



(21) 



dr 



Substituting this expression into ([20)1 yields the follow- 
ing solution for the disk radiative flux: 



8^ 
3 /riN2 



1 



Ma(r) 



4 \ r 

m 

87rr2 



dMa 



47rr 



(22) 



dr 



The first two terms on the right hand side of this equa- 
tion describe the rate at which gravitational binding 
energy is extracted from matter as it accretes (i.e. loses 
angular momentum); the second term in particular de- 
scribes the rate at which nonzero MHD stresses at the 
innermost stable orbit locally dissipate turbulent en- 
ergy (in practice, a convenient inner boundary condi- 
tion is to set this term equal to zero, although in prin- 
ciple, energy can still be extracted beyond this bound- 
ary). The last term on the right hand side of de- 
scribes the rate at which energy is removed from the 
disk by outfiows involving mass and Poynting fluxes. 

The result ([22)l for the radiative flux of an accre- 
tion disk corrected for the effects of outflows can be 
expressed as 



where 



3GMMa(r-) 
87rr^ 



(19) Mr) = 



1 - 



Ma(n) 

Ma(0 



n\i/2 

r 



(23) 



(24) 
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is the accretion factor and 



/w(r) 



1 



Ma(r)r2ri Jr 



A 2/, \- 

w^r— J- 47rr {t^^) 



dr 



(25) 



is the outflow correction factor (the 'w' subscript de- 
notes wind). Note that /w is equivalent to the frac- 
tional rate of vertical angular momentum transport in 
the disk. From eqns. (fT5|) and (fT6|) . we have 



/w(0 = ^ 

Ur) 
Similarly, 



Ur) = 1 



and hence. 



ja(n) + ir(n) 



(26) 



(27) 



/a(r)-/w(r) 



ia(r) 



Therefore, the efficacy of disk radiant emission is en- 
tirely determined by the relative importance of radial 
transport of angular momentum in the disk. If vertical 
transport of angular momentum dominates, then en- 
ergy is efficiently removed from the disk before being 
locally dissipated therein. 

It is noteworthy to compare our outflow-modified 
disk flux with the disk flux predicted by standard ac- 
cretion disk theory. In the absence of outflows (that 
is, when dMa/dr = and {t^z)~^ = 0), then = 0, 
Ma(ri) = Ma(r) = M and hence, /a = 1 - (n/U^^. In 
the limit of vanishing stresses at ri, this is equivalent 
to the small— r c orrection factor from stand ard accre- 
tion disk theory (jShakura fc SunvaevI 119731) . In other 



words, standard accretion disk theory predicts that an- 
gular momentum is transported radially outwards only 
and consequently, all the extracted gravitational bind- 
ing energy is locally dissipated in the form of radiation. 



3 Disk spectrum 

An expression for the total disk radiative luminosity, 
Ld, is obtained by integrating each term in eqn. ((^^ 
over all disk radii, from r = Vi to r = oo. Similarly, 
the spectrum of disk emission can be calculated assum- 
ing local blackbody emission and summing the spec- 
trum from each annulus. This is the Multi-Colour Disk 
(MCD) model used for standard disks (|Mitsuda et al.l 
Il984l) . We will henceforth refer to our generalized disk 



prescription as the Outflow-Modifled Multi-Colour Disk 
(OMMCD) model. 

Before we can calculate the total luminosity and 
emission spectrum for the OMMCD model, it is nec- 
essary to define a specific functional form for the out- 
flows. Realistically, of course, whether a disk becomes 
conducive to outflows and the resulting behaviour of 
those outflows are determined by several factors, in- 
cluding the degree of fleld-matter coupling, as well as 
the field strength and topology, all of which will in gen- 
eral vary with radius across the disk. Indeed, numerical 
simulations to date have been unable to show whether a 
nonzero, nett poloidal field can evolve from a stochastic 
field with zero nett flux. Here, we wish to obtain a sim- 
ple quantitative measure of how outflows might affect 
the observable properties of an accretion disk. 

3.1 A Simple Model 

A simple phenomenol ogical model for the ni ass accre- 
tion rate is (see e.g 



Wardle fc K5^ll99i H llQQd 



(28) Casse fc Ferreira 200C for other self-similar models) 



Ma(r-) 




r < r« 
r > 



(29) 



where represents a critical wind radius beyond which 
the mass outflow is negligible. The wind mass loss rate 
then satisfles Af„(r) = M — Ma(r), from the continuity 
equation (fTT|) . We can now use this model to explic- 
itly calculate the factor /a, which appears as a source 
term in the disk flux solution (j23|) and which is deflned 
by ((HI). For the simplest model, we assume that the 
stresses Tr^{ri) at the inner boundary radius vanish, so 
that 




r < r„ 
r > r„ 



(30) 



Similarly, we can define a specific model for fw{r), 
given by (|25p . From the vertically-integrated, differen- 
tial form of the angular momentum equation, p3[) . we 
have 



.dMa 
dr 



d 



(31) 



-f 

dr V 



i:iM.r^^) 



We wish to define a model in which the relative im- 
portance of the outflows in transporting angular mo- 
mentum decreases with increasing radius. We choose a 
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simple power-law decline, such that 



(32) 



with w > 0. Thus, in the limit w —>^ 0, radial transport 
of angular momentum by the internal MHD stresses 
Tr^ is negligible and the dominant mode of transport is 
via outflows. Inserting the relation ([5^ into PT|) and 
simplifying, we have 



dr 



-(1/2 



This implies the following relation for the outflow factor 
/w defined by ([25|l: 



i/2+p-wJ Vn 

1- ( - 



/w(?' < r-w) = 



(34) 



fwir > rw) 



1/2+p 



1/2 + p — w / \r. 

(1/2+p— uj) 

1 - 



N -1/2 , 

7'i 



1/2 



1/2 



l-(^ 



r 

w / \ri 

-ii/2-wy 



(35) 



Eqns. ((5U]l and for /a and /w, respectively, now 
specify the solution for the OMMCD radiative flux as 
a function of radius, given by eqn. (|22p . 

The total disk luminosity can now be calculated by 
integrating the flux over all disk annuli, as follows: 



where 



1 GMM (l + 2p 



1 



3p 



(1 + 2p) 



p 



(38) 



is the total accretion power and 

l + 2p 



_ 1 GMM 



1-p 



-p 



(33) 

-p+w) 



p(3 + 2w) 
{l + 2p){l + w) Vn 



(39) 



is the total power removed by outflows. Note that in 
the limit p — » 0, corresponding to negligible mass loss 
(i.e. Ma —^ constant), the accretion power approaches 
the solution from standard accretion disk theory. Pa = 
GMMa/2ri, and the wind power is Pw — * ^a/(l + w). 
Thus, a substantial fraction of the total accretion power 
can be removed from the disk by a magnetic torque 
alone. This solution represents a Poynting flux domi- 
nated outflow and can be identified with the formation 
of relativistic jets that carry away very little matter, 
but transport a large amount of mechanical energy very 
efficiently across vast distances. 

The emission spectrum predicted by the OMMCD 
model can be calculated in the same way as that for 
an MCD model (i.e. a standard accretion disk). As- 
suming each annulus in the disk radiates locally like a 
blackbody, P^, the total disk spectrum is calculated by 
summing the individual spectra from each annulus: 



Ld = 2 / Fd{r)2Trrdr 
3 GMM f°° Ma(r) 



2 n Jr, M 
We can write the solution as 



(37) 



Ld,v — 2 



7rP^[T(r)]27rrdr 



(40) 



(36) 



[/a(r)-/w(r)] - 

n 



where T{r) = [Fd{r)/aY/'^ is the efi^ective disk temper- 
ature of each annulus and a is the Stefan-Boltzmann 
constant. 

Figure [T] shows theoretical disk emission spectra 
Pd.iv for three different models of black hole disk ac- 
cretion: 1. the MCD model (blue dotted curve) de- 
rived from the standard "Shakura-Sunyaev" theory 
(jShakura fc SunvaevI [ToTsI) : 2. the OMMCD model 
with negligible mass loss in an MHD outflow (p = 0.1, 
w — 0.1, red dashed curve); and 3. an OMMCD model 
with stronger mass loss in an MHD outflow {p — 0.5, 
w = 0.1, red solid curve). In both OMMCD models, 
the mechanical energy carried away by the outflow is 
« 90 % of the total accretion power and mass loss is 
restricted to radii r < lOOrg. In all three models, M = 
IQ^ Mq (corresponding to pEdd ~ 1-3 x lO^^ergs^^) 
and M — O.27M0yr~i. This is the nett mass influx 
at r = oo and corresponds to the Eddington rate when 
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14.0 14.5 15.0 15.5 16.0 16.5 17.0 

log [V (Hz)] 

Fig. 1. — Comparison of theoretical accretion disk spec- 
tra . The MCD model (dotted blue curve) is based on 
standard accretion disk theory. The OMMCD model 
(red curves) predicts a radial disk structure modified 
by outflows. The dashed red curve is the predicted 
OMMCD spectrum for a Poynting flux dominated out- 
flow, while the solid red curve is the OMMCD spectrum 
for a mass-loaded outflow. In all models, the black hole 
mass is M = lO^M© and the mass accretion rate at 
infinity is M = 0.27 Mgyr^^ 

the radiative efficiency is 77 ~ 0.08. This is strictly only 
the case for the MCD disk model. When outflows are 
important, the radiative efficiency can fall below the 
nominal 8% predicted by standard theory because the 
outflows remove energy from the disk (as well as an- 
gular momentum), so there is relatively less energy to 
dissipate and radiate away. For the OMMCD models 
shown in Fig.[Tl the disk radiative efficiency is about an 
order of magnitude lower than that of a standard disk. 

As is clearly evident in Fig. [I] outflows can substan- 
tially modify the overall spectrum of emission. The gen- 
eralized OMMCD model predicts a high-energy cut-off 
in the emission spectrum at lower energies than that 
predicted by the standard MCD model. This is a di- 
rect results of the importance of outflows at small radii, 
where the highest-energy emission originates. In the 
case of stellar-mass black holes, the spectrum is most 
affected at X-ray energies (around a few keV), while 
for supermassive black holes, it is the extreme ultravi- 
olet region 10^^ Hz) that is most affected by out- 
flows. The OMMCD model also predicts a broadband 
region of the disk spectrum (typically around visual 
wavelengths) that is much flatter than the character- 
istic v^/^ law predicted by standard theory. 



We have used a specific phenomenological model to 
quantitatively calculate how a disk spectrum can be 
modified by outflows and a different choice for the func- 
tional form of the outflows may result in a somewhat 
different predicted spectral shape. Nevertheless, the 
overall effect of outflows on a disk should be largely in- 
sensitive to the detailed physical properties of the out- 
flows; the extraction of angular momentum and binding 
energy from the disk will inevitably result in an emis- 
sion spectrum that is dimmer and redder than that pre- 
dicted by the standard model for the same mass accre- 
tion rate. 

3.2 Comparison with observations 

Our results demonstrate that the modifications to a 
standard disk spectrum as a result of energetic outflow 
phenomena are clearly not negligible when the outflows 
are primarily responsible for the removal of angular mo- 
mentum from accreting matter. Direct observational 
verification of the predicted OMMCD spectrum may 
be possible for only some types of accreting sources, 
however, as the outflows will also produce their own 
characteristic radiative signatures that may overlap in 
the spectral energy band where we predict the disk to 
be most strongly modified. This can confuse our inter- 
pretation of the emergent observed spectrum. 

In the case of galactic X-ray binaries (XRBs), 
for instance, the disk spectrum can be substantially 
modified at X-ray energies (typically around 1 keV) . 
These sources often exhibit a power-law X-ray spec- 
trum above 1 keV that is thought to arise in a mag- 
netized, tenuous corona and/or relativistic jet as a re- 
sult of inverse Compton scatter ing off disk photons (see 
McClintock fc RemillardI I2OO6I for a review). Existing 
models for this emission process simply use a standard 
disk spectrum for the seed photon distribution. How- 
ever, this is clearly not a self-consistent calculation since 
the formation of a corona and/or jet from disk magnetic 
fields inevitably modifies the disk radial structure and 
hence, photon spectrum, as we have shown here. Thus, 
although we would not expect to directly see the modi- 
fied disk spectrum in XRBs, our model can be indirectly 
tested by using the underlying OMMCD spectrum as 
the source of seed photons that are upscattered in a 
disk corona and/or jet. By comparing the predicted 
X-ray spectrum against the observed spectrum, we can 
improve our interpretation of the source characteristics 
and obtain tighter constraints on key physical parame- 
ters such as mass accretion rate. 

The above type of calculation has been p erformed by 
us for ultra-luminous X-ray sources (ULXs) (jFreeland et al.^ 
2006[ ). These exceptionally luminous XRBs probably 
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involve accretion onto a black hole more massive than 
those found in galactic XRBs. They provide an ef- 
fective test for our modified disk model because some 
ULXs exhibit a low energy ( < 1 keV) spectral compo- 
nent that can be interpreted as disk emission. How- 
ever, because these sources emit virtually all of their 
radiative power as hard, power-law X-rays, which must 
necessarily be produced outside the disk, the accretion 
disk must be substantially modified and a standard 
disk model should not be use d to fit the soft sp ectral 
component. Indeed, we show ( Kuncic et al. 2006f ) that 
the OMMCD spectrum can adequately fit the observed 
soft spectral component, implying a black hole mass 

M - lOOAfg. 

A similar test can also be performed for quasars and 
other AGN, powered by accretion onto a supermassive 
(106-9 Mq) black hole. The prominent optical/UV 
continuum feature known as the "big blue bump" 
( Sanders et al. 19891 ) is generally interpreted as ac- 
cretion disk emission (IShieldslll97i iMalkan fc SargentJ 
19821) . However, attempts at fitting accretion disk spec- 
tral models to the observed spectr a have had mixed 
success (see Koratkar fc Blaes 19991 for a review). For 
AGN, the disk spectrum will be most strongly mod- 
ified by outflows at UV wavelengths (see Fig. [T]). 
Although there is less overlap between the modified 
disk and corona/jet X-ray emission than in the XRB 
case, the poor transmission of UV radiation through 
Earth's atmosphere means that direct observations 
of the predicted OMMCD spectrum in AGN are not 
straightforward. Nevertheless, spectral energy distri- 
butions (SEDs) have been compiled fro m observations 



of high-redshift quasar s with the SDSS ((Richards et al 



20061 iTrammell et al.l 120071) and from satellite UV 



Telfer et al.l 120021) and with FUSE 



Shang et al.l I2OO5I ) 



Zhene et al. 


1997 


(IScott et al. 


2004 



Interestingly, these observations 
hav e revealed a far-UV b reak in the SEDs of quasars 
(see ^Trammell et al. 2007 and references cited therein) . 
Specifically, the observed spectra decline dramatically 
at wavelengths shorter than 1100 A (i.e. frequencies 
above 10^^'^ Hz). This is difficult to explain with stan- 
dard accretion disk models, which predict that the spec- 
trum should continue to rise into the far and extreme 
UV (see Fig. p. 

Figure [2] shows a comparison of the MCD and 
OMMCD spectra against the in ean SED of SDSS type 
I quasars ( Richards et al.|[2006l ). The parameters used 
for the MCD spectrum (dotted blue curve) are: M = 
1 X IO^Mq and M = = 2.7Mqyy-^, which is 
equivalent to O.lMEdd- The OMMCD spectrum uses 
M = 5 X IO^Mq and M = 25M0yr-\ which corre- 
sponds to a mass accretion rate of about 19MEdd at 
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Fig. 2. — A comparison between observed and pre- 
dicted quasar optical/UV luminosity spectra. The solid 
black curve is the restframe mean spectral energy distri- 
bution of SDSS type I quasars, corrected for galactic ex- 
tinction (see Richards et al. 2006 for online data) . The 
dotted blue curve is the spectrum predicted by a stan- 
dard multi-colour disk (MCD) model with M = IO^Mq 
and Ma = 2.1 M^yx-^ (equivalent to = O.lMEdd)- 
The dashed red curve is the spectrum predicted by the 
outflow modified multi-colour disk (OMMCD) model, 
with M = 5 X IO^Mq and M = 25MQyr-i (corre- 
sponding to A/a(oo) = 19MEdd), where 90% of the ac- 
cretion power is carried away by an MHD outflow. 

r = 00. The outflow parameters are p = 0.5, w = 0.1 
and Tw/n = 10. In this case, only 30% of the matter 
accreting at r = cxd reaches r\ and 90% of the total 
accretion power is removed by the outflow. Because 
most of this energy is removed preferentially at small 
radii, the disk is not as hot at these radii as a stan- 
dard disk. Consequently, the emergent disk spectrum 
is fainter in the far-UV, precisely in the bandpass where 
a break has been detected independently by several in- 
struments. Note, however, that both the MCD and 
OMMCD models underpredict what appears to be a 
secondary spectral component in the extreme UV (at 
frequencies above 10^^'^ Hz). This component resem- 
bles a power-law with a spectral index close to 2 and 
could result from Comptonization. 

Fig. [2] clearly shows the difference between the far- 
to-extreme UV spectrum predicted by a standard disk 
and that predicted by an outflow modified disk. There 
are two important points to draw from this compar- 
ison. Firstly, the standard model does not take into 
account the partitioning of accretion energy into disk 
radiation and other non-radiative forms of energy (e.g. 
Poynting flux) . The second important point is that the 
two different disk models predict different values for the 
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fundamental physical parameters M and M. The impli- 
cation is then that we could be severely misinterpreting 
quasars and AGN by using the standard disk model to 
infer, in particular, mass accretion rates. The univer- 
sality of accretion related outflow phenomena suggests 
that the true rates of mass accretion could be consid- 
erably higher than we are currently estimating them to 
be. 



4 Concluding Remarks 

The standard theory for astrophysical disk accretion 
has enjoyed remarkable success since its inception over 
30 years ago. However, ongoing rapid advances in in- 
strumental technology are placing increasingly tighter 
observational constraints on interpretational models 
that utilize the standard theory. The ubiquity of en- 
ergetic outflow phenomena in accreting astrophysical 
sources gives us two important clues to the underly- 
ing processes responsible for accretion: (1.) Vertical 
transport of angular momentum cannot be ignored and 
may indeed be the dominant mode of transport at small 
radii; and (2.) The removal of accretion energy from the 
disk results in a modified radial disk structure at small 
radii and a disk spectrum that can be substantially dif- 
ferent from that predicted by standard accretion disk 
theory. 

We have presented a simple model that can be used 
to quantitatively calculate model spectra for disks mod- 
ified by a magnetized jet and/or corona and by mass- 
loaded winds. These model spectra can be applied to a 
variety of different sources to improve our interpetation 
of the observations and more accurately determine key 
physical parameters such as mass accretion rate and 
black hole mass, as well as the partitioning of accretion 
power into radiative and non-radiative forms. 

The authors wish to thank all the particpants for 
contributing to a very successful and very informative 
Fifth Stromlo Symposium. 
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